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摘 要 : 对 于 复杂 的 体系 ， 组 成 体系 的 各 粒子 间 的 相互 作用 决定 其 能 量 状态 。 对 于 不 同 势 场 ， 描 述 粒 子 运动 规律 的 Schrodinger 
方程 是 不 同 的 微分 方程 。 本 文 在 原来 的 Hamilton 谐振 子 数 学 模型 基础 上 ， 通 过 引入 Morse 势能 函数 代替 原来 的 简 谐 振子 势 模 
型 ， 研 究 了 一 维 气 体 分 子 链 中 声 传播 的 问题 。 从 能 量 的 概念 出 发 ， 将 声 传播 问题 转换 成 求解 谐振 子 的 波 函 数 以 及 能 量 本 征 值 的 
问题 ， 建 立 配 分 函数 与 热力 学 量 之 间 的 关系 ， 以 Morse 势能 函数 的 振子 模型 构造 的 声 传播 模型 求 取 声 压 及 声 能 量 。 通 过 比较 ， 
经 典 方法 与 量子 谐振 子 模型 计算 的 一 维 声 压 吻合 。 
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Sound Propagation Research on the Oscillator Model of Morse Potential 


Energy Function 
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Abstract: For a complex system, the reciprocal function of the particles constituting the system decides its energy 
state. Schrodinger equation which describes the particles moving rule will vary im different potential field. In this 
paper based on the mathematic model in Hamilton harmonic oscillator researched before, Morse potential energy 
function which replace the harmonic oscillator potential model is introduced to study sound propagation in the 
one-dimensional gas molecule chain. From the standpoint of energy, the problem about sound propagation is 
converted into the problem of the harmonic wavelet function and the energy Eigen values. The relationship between 
the partition function and the thermodynamics quantity is built and sound propagation model based on oscillator 
model of Morse potential energy functions is used to get sound pressure and sound energy. Finally the 
one-dimensional sound pressure results got by quantum oscillator model agree well with classical approach results. 
Keywords: aeroacoustics; Morse potential energy function; Schrodinger equation 


0 引言 


已 有 气动 声学 理论 研究 基本 上 是 在 牛顿 体系 1 李 政道 提 到 门 :“20 世纪 的 文明 是 微观 的 ，…… 
展开 研究 的 ， 虽 在 工程 应 用 上 获得 巨大 成 功 ， 但 它 ”我 认为 21 世纪 微观 和 宏观 应 结合 成 一 体 。……… 昌 
们 从 建立 之 初 就 有 其 不 足 之 处 ， 令 其 发 展 缓慢 ， 近 ”前 微观 和 宏观 的 冲突 己 经 非常 尖锐 ， 靠 一 个 不 能 解 
几 十 年 其 或 停滞 不 前 、 难 以 改善 。 人 们 不 仅 要 在 目 。 决 男 一 个 ， 把 它们 联系 起 来 会 有 一 些 突破 ， 这 个 突 
前 的 手段 开展 气动 噪声 研究 ， 而 且 还 要 考虑 在 未 来 。” 破 会 影响 科学 的 未 来 。” 如 有 一 个 公共 的 力学 体系 
的 计算 条 件 下 开展 气动 声学 问题 的 研究 。 目 前 无 论 。 能 将 气动 声学 的 微观 和 宏观 研究 联系 起 来 ， 为 气动 
是 工程 应 用 还 是 理论 研究 ， 要 求 气动 声学 理论 能 提 ” ”声学 研究 找 出 一 种 新 途径 ，Hamilton 力学 就 有 这 样 
供 声 波 运动 及 声 流 相 互 作 用 规律 ， 以 及 在 理论 上 提 优点 。 Hamilton 力学 就 具有 这 样 优点 。 如 有 一 个 公 
供 了 更 为 清晰 的 物理 图 象 。 共 的 力学 体系 能 将 气动 声学 的 微观 和 宏观 研究 联系 
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起 来 ， 为 气动 声学 研究 找 出 一 种 新 途径 ，Hamilton 
力学 就 有 这 样 优点 。 从 某 种 意义 上 讲 ， 哈 密 顿 体系 
是 更 本 质 和 更 重要 的 ， 一 切 保守 体系 都 可 以 在 哈密 
顿 体系 下 统一 表示 。 
文献 中 论述 了 在 Hamilton 系统 下 ， 以 谐振 子 势 
能 模型 结合 群 论 建立 数学 模型 讨论 了 一 维 气体 分 子 
0 题 。 

对 于 复杂 的 体系 ， 组 成 体系 的 各 粒子 间 的 相互 
作用 决定 其 能 量 状态 。 各 相互 作用 以 及 处 理 中 中 心 
场 使 用 近似 方法 ， 使 得 复杂 体系 的 Schr5dinger 方程 
中 有 可 能 出 现 高 次 非 谐振 子 势 、 极 化 等 效 势 、 环 形 
非 球 振子 势 、 势 等 高 次 正 蝴 与 逆 蝴 势 函 数 以 及 它们 
的 有 登 加 等 不 同 的 势 场 。 对 于 不 同 势 场 ， 描 述 粒 子 运 
动 规律 的 Schr6dinger 方程 是 不 同 的 微分 方程 ， 但 只 
有 少数 势 场 的 Schr6dinger 方程 才能 得 到 其 特征 值 的 
解析 解 ， 多 数量 子 体系 的 Hamilton 算数 比较 复杂 ， 
Schr5dinger 方程 不 能 严格 求解 品 。 
任何 体系 在 平衡 位 置 附近 的 小 振动 ， 例 如 分 子 


振动 、 晶 格 振 动 、 原 子 核 表 面 振动 以 及 辐射 场 的 
动 等 往往 都 可 以 分 解 成 若干 彼此 独立 的 一 维 简 谐 
动 中 。 简 谐振 动 往往 还 作为 复杂 运动 的 初步 近似 四。 
对 于 双 原 子 分 子 谐振 而 言 ， 应 用 上 述 简 谐振 子 
模型 可 以 较 好 地 近似 描述 双 原 子 分 子 在 低能 量 时 的 
谐振 行为 。 但 是 在 实际 上 ， 由 于 双 原 子 分 子 振动 并 
不 是 简 谐 振动 ， 振 动能 级 也 并 不 均匀 ， 在 常温 下 双 
原子 分 子 几 乎 处 于 低能 级 上 ， 这 时 双 原 子 分 子 势能 
曲线 与 简 谐 振子 的 势能 曲线 在 低能 级 振动 时 相近 门 。 
本 文 在 原来 的 Hamilton 谐振 子 数学 模型 基础 
上 ， 通 过 引入 Morse 势能 函数 代替 原来 的 简 谐 振子 
势 模 型 ， 研 究 了 一 维 气体 分 子 链 中 声 传播 的 问题 。 
从 能 量 的 概念 出 发 ， 将 声 传播 问题 转换 成 求解 谐振 
子 的 波 函 数 以 及 能 量 本 征 值 的 问题 外 。 本 文 与 文献 站 
的 不 同 在 于 : 分 子 间 只 与 相 邻 的 分 子 相互 作用 ， 取 
Morse 势能 函数 , 数值 求解 分 子 控制 方程 时 用 四 阶 辛 
算法 ， 同 时 比较 了 经 典 方法 与 量子 谐振 子 模型 计算 
的 一 维 声 压 。 


1 模型 建立 的 理论 基础 


本 文选 取 一 维 体系 进行 研究 ， 空 气 成 分 中 占 绝 
大 部 分 的 是 氮气 和 氧气 ,他们 都 是 双 原 子 分 子 中 。 分 
子 间 互 作用 原理 详 见 文献 哲 。 

可 证 明 长 波 极限 下 ， 一 维 分 子 链 的 振动 可 化 为 
连续 体 弹 性 波 。 一 维 分 子 链 的 运动 方程 为 : 


展 
展 
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m 谎 = LOL 十 Hn 和 21,) 
[La 一 他 Hn Hn ( 1 》 


= Ba’ a a 
a 


在 长 波 极限 下 ， 位 移 可 以 看 作 连 续 变 化 的 ， 即 : 
1 区 = Pa’ (2) 
而 弹性 波 方 程 为 : 
B= vi, (3) 
对 比 得 波 速 为 : 
v= je (4) 
m 
设 试探 解 为: 
u, = Ae ("4) (5) 


式 中 4 为 振幅 ，w 为 圆 频率 ，g 为 波 矢 ， 与 波长 
4 的 关系 为 : q=2x/Anag 为 第 n 个 分 子 振动 的 位 相 ， 
将 〈2) 代入 (1)， 容 易 求 得 w 与 9 的 关系 为 : 
,_280 


0 二 


—cosga) (6) 


由 式 (6) 不 难看 出 ， 当 na-ma= 信 ， 即 第 n 和 第 
m 个 分 子 的 位 移 相 等 ， 所 以 式 〈6) 所 描述 的 原子 转 
绕 平 衡 位 置 运动 是 以 行 波 的 形式 在 气体 中 传播 的 ， 
是 气体 介质 中 分 子 子 的 一 种 集体 运动 形式 。 
色散 关系 可 得 到 波 的 相 速 度 和 群 速度 为 : 
相 速 度 : 


. Ga 
$ sin — 
y -= 2 (7) 


q m) gaa 
2 
群 速度 : 
疙 
“= 和 -| 人 i (8) 
dq m 2 
由 此 可 以 看 到 ， 由 于 分 子 的 不 连续 性 ， 波 的 相 
a 
数 ， 这 是 因为 当 波长 很 大 时 ， 一 个 波长 范围 含 者 干 


分 子 ， 相 邻 分 子 的 位 相差 很 小 ， 分 子 的 不 连续 效应 
很 小 ， 波 接近 于 连续 媒质 中 的 弹性 波 。 


2 计算 模型 


文献 叫 论 述 了 在 Hamilton 系统 下 ， 以 谐振 子 势 
能 模型 结合 群 论 建立 数学 模型 讨论 了 一 维 气体 分 子 
链 中 声 传播 的 问题 。 本 文 基于 文献 中 基础 上 , 程序 编 
译 思想 如 下 图 1 所 示 。 
基于 群 论 及 谐振 子 模型 的 声 传播 方法 步 又 为 : 

1、 先 将 声波 传播 问题 转换 成 求解 Schr5dinger 
方程 ， 确 定 波 函 数 及 能 量 本 征 值 〈 能 级 )。 这 里 用 谐 
振子 的 Schrodinger 方程 ， 后 续 研 究 可 以 再 考虑 复杂 
体系 县 加 势 的 含 时 Schrodinger 方程 ; 

2、 求 解 分 子 正 则 运动 方程 ， 解 久 期 方程 ， 得 振 
动 频率 及 振幅 ; 

3、 通 过 声 能 与 声 压 的 关系 获得 声场 有 关 量 。 


图 1 程序 编译 思想 图 
Fig.1 Chart of program compile though 
3 分 子 Morse 势 能 函数 谐振 子 配 分 函数 
对 于 双 原 子 分 子 谐振 而 言 ， 应 用 简 谐 振子 模型 
可 以 较 好 地 近似 描述 双 原 子 分子 在 低能 量 时 的 谐振 
行为 . 但 是 在 实际 上 , 由 于 双 原 子 分 子 振动 并 不 是 简 
谐振 动 ， 振 动能 级 也 并 不 均匀 ， 在 常温 下 双 原 子 分 
子 几乎 处 于 低能 级 上 ， 这 时 双 原 子 分 子 势能 曲线 与 
简 谐 振子 的 势能 曲线 在 低能 级 振动 时 才 相近 。 图 2 
表示 谐振 子 、Morse 势能 函数 振 势 能 曲线 ， 图 3 表示 
谐振 子 、 Morse 势 振子 分 子 间作 用 。 
这 里 用 Morse 势能 函数 方法 ， 它 可 以 更 好 地 接 
近 双 原子 分 子 谐振 的 实际 情况 .下 面 对 此 加 以 讨论 . 
引入 双 原 子 分 子 势能 函数 


V(r)=D,|1-exp(-a(r—x,))] (9) 


式 中 : 


D, 一 一 分 子 从 势能 曲线 最 低 点 到 孤立 原子 


(r —> %) 时 的 能 量 ; 


六 一 一 双 原 子 的 平衡 核 间 的 距离 ; 
4 一 一 与 分 子 有 关 的 参量 。 
当 (x 一) 很 小 时 略 去 其 高 阶 项 ， 就 可 以 利用 势能 


函数 级 数 展开 式 求 解 Schr5dinger 方程 。 这 里 ， 需 要 
确定 能 量 零点 ， 可 分 别 取 势 能 曲线 最 低 点 或 者 分 子 
振动 基态 为 能 量 零点 作为 能 量 零 点 。 


=H 


对 势能 函数 V(r) 平衡 点 做 Taylor 级 数 展开 为 


PO)=7(n)+[ 竺 ] CC- 
1{ ar 2 (10) 
导入 9 
1 1 d 
站 | (r—z.,) +K 
如 果 取 势能 曲线 最 低 点 作为 能 量 零点 ， 可 以 得 


到 它 的 振动 能 级 及 其 相应 的 配 分 函数 分 别 为 : 


1 be 
E=| n+t+~ hv,—lnt+e | xhv, 
2 2 
之 
万 h 
式 中 v= 上 x, Vo 
x \24 4D, 


可 由 光谱 数据 求 得 。 
如 果 取 分 子 振动 基态 为 能 量 零点 ， 则 振动 能 级 
与 相应 的 配 分 函数 可 写成 如 式 〈12) 所 示 。 

已 =| zza(a+lDz |hv, 


， hy (12) 
g' = Zeo|-D-o0 +1)x, 他 | 
由 于 双 原 子 Morse 势能 函数 中 的 D, 是 双 原 子 从 
平衡 位 置 离 解 成 2 个 相距 无 限 远 的 原子 所 需 的 能 量 ， 
此 振动 能 级 已 , 最 大 为 万 ,与 此 相应 的 求 和 的 分 子 


EH 


数 n 应 取 为 上 限 n,。、， 它 由 式 (13) 确定 ， 即 


max 


2 
万 = 区 3m 二 hn 
2 2 


将 ze 和 代入 可 得 1 1 


(13) 


q' = 六 exp -[n —n(n +D 全 | 
n=0 


201802.00070v1 


InaXIV: 


ch 


式 中 : 0， 双 原 子 分 子 振动 特征 温度 ， 它 等 于 
六 (15) 
Kk 


2 谐振 子 、 Morse 势能 函数 振 势 能 曲线 
Fig.2 Vibrational potential energy curve of harmonic 


oscillator and Morse potential energy function 
Us T 


Ba- 7 入 


:CA EC ZSI JI ZE 
3 谐振 子 、Morse 势 振子 分 子 间作 用 
Fig.3 Interaction of harmonic oscillator and Morse potential 


Vibrator of molecule 
4 能 级 与 宏观 量 的 联系 
如 果 已 知 气体 的 自由 能 函数 (7T,V)， 六 为 体 
积 ， 就 可 以 根据 : 
p=-(); Cle 


OV 
写 出 气体 的 状态 方程 。 自 由 能 函数 可 以 一 般 地 


F=-kTInZ (17) 
2 为 配 分 函数 : 
-E, 
a hr (18) 


连 加 式 是 对 所 有 分 子 的 能 级 , 相 加 。 


能 级 忆 , 除 包括 分 子 处 于 格 点 时 的 平衡 气体 的 能 


量 U(V) 外 ,还 有 各 波 的 振动 能 于 Cu + ho， 7 


J 
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标志 各 不 同 波 ，n; 为 相应 的 量子 数 。 配 分 函数 Z 包 


muy 


系统 的 所 
2，… 相 加 ， 从 而 得 到 ; 


pa eV Je I bb 0 


/ 


二 or 1 
—ho, 
le hr 


代入 自由 能 公式 (9) 得 到 : 


量子 态 ， 因此 应 分 别 对 每 个 n, =0， 1， 


ni=0 


(19) 


FUtry me (20) 
当 体系 体积 改变 时 , 波 频率 也 将 改变 , 所 以 上 式 除 U 
以 外 ,各 频率 w 也 是 宏观 参量 的 函数 。 根据 (12) 
对 信 求 微 商 ， 得 到 : 

dU 1 n |]a 
[es ?3 此 (21) 


J KT-1 


上 式 包含 了 各 振动 频率 对 的 矿 依赖 关系 ， 因 此 具有 
很 复杂 的 性 质 。 格 律 乃 森 针 对 这 种 情形 ， 提 出 一 个 
有 用 的 近似 。 如 把 上 式 写 成 : 


a Tl2 fa) dr 


其 中 式 (22) 括号 内 是 平均 振动 能 。 


dlnow., 
式 Q2) 中 表征 频率 随 体积 变化 的 一 一 也 是 一 个 
n 


无 量 纲 的 量 , 假设 它 近 似 对 所 有 振动 相同 , 这 样式 (22) 
就 简化 为 下 列 近似 状态 方程 : 


站 = 一 -一 +7Y 二 (23) 


革 中 万 表示 气体 的 平均 振动 能 : 


0 (24) 


dlnV 
7 为 常数 。 一 般 @ 随 增 广 加 而 减 小 ，y 上 共有 正 
的 数值 。 


5 声 传 播 模 拟 及 分 析 
空气 中 15'C 时 声速 为 340 m/s， 假 设 各 个 振子 


振 前 都 处 于 基态 ，n, = 0 。 这 里 ， 由 于 同 核 双 原子 分 


子 N, 、O0, 属 于 Dw 群 ,计算 时 以 D。 群 操 作 。 简 便 


起 见 ， 选 取 点 震源 为 <= sin(0.002zn0 ， 则 所 考察 的 系 
统 的 总 能 量 为 : 
uM 1 N 
B= ho n+) = ho (25) 


从 能 量 表达 式 看 出 ， 各 个 振子 的 波 数 或 频率 可 
以 不 一 致 ， 故 一 系列 谐振 子 的 组 合 相当 于 波 包 。 

模拟 过 程 中 把 波动 过 程 看 成 是 一 种 能 量 的 流 
动 。 即 声波 对 谐振 子 做 功 ， 体 系 就 会 发 热膨胀 。 图 4 
为 谐振 子 模型 模拟 声波 传播 过 程 中 ， 声 压 与 密度 变 
化 关系 图 。 图 5 表示 的 是 用 二 阶 有 限 差分 法 FDTD) 
模拟 的 经 典 声 传播 过 程 中 ， 声 压 随 波 程 (时 间 与 速 
度 的 乘积 ， 波 速 为 常量 ) 变化 示意 图 。 通 过 比较 ， 
新 模型 声 压 随 波 程 变化 的 趋势 和 经 典 声 传播 的 趋势 
十 分 吻合 。 这 说 明 本 文 数值 计算 有 效 。 


压 和 密度 变化 


波 程 x 


图 4 声 压 与 密度 随 波 程 变化 比较 图 
Fig.4Comparison chart between sound pressure and density 
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图 5 一 维 二 阶 FDTD 经 典 声 传播 声 压 示 意图 
Fig.5 Diagram of one-dimensional second order FDTD classic 


sound propagation sound pressure 


6 结论 


本 文 基于 用 Morse 势能 函数 的 振子 模型 及 群 论 
方法 对 声 在 空气 中 传播 建立 了 一 个 新 的 哈密 顿 数 学 
模型 。 分 析 并 讨论 分 子 状态 下 ， 分 子 能 级 与 配 分 函 
数 的 关联 ， 并 通过 配 分 函数 与 热力 学 量 之 间 的 关系 ， 
以 Morse 势能 函数 的 振子 模型 构造 的 声 传播 模型 求 
取 声 压 及 声 能 量 。 通 过 比较 ， 经 典 方法 与 量子 谐振 
子 模型 计算 的 一 维 声 压 吻合 。 

此 外 ， 新 模型 只 考虑 了 声波 传播 问题 ， 将 问题 
转换 成 求解 谐振 子 的 波 函 数 以 及 能 量 本 征 值 的 问 
题 ， 除 了 离散 关系 和 本 征 值 方程 的 代数 操作 , 它 既 不 
需要 完整 求解 传播 模型 ,也 没有 任何 迭代 计算 。 通 过 
与 有 限 差分 法 (FDTD) 模拟 结果 的 比较 ， 新 模型 的 
数值 计算 有 效 。 
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